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ABSTRACT 

(N ; 

Context. The WASP-10 planetary system is intriguing because different values of radius have been reported for its transiting exoplanet. 
' The host star exhibits activity in terms of photometric variability, which is caused by the rotational modulation of the spots. Moreover, 

X^/^ ' a periodic modulation has been discovered in transit timing of WASP-10 b, which could be a sign of an additional body perturbing 

the orbital motion of the transiting planet. 

Aims. We attempt to refine the physical parameters of the system, in particular the planetary radius, which is crucial for studying the 
internal structure of the transiting planet. We also determine new mid-transit times to confirm or refute observed anomalies in transit 
timing. 

Methods. We acquired high-precision light curves for four transits of WASP-10 b in 2010. Assuming various limb-darkening laws, 
we generated best-fit models and redetermined parameters of the system. The prayer-bead method and Monte Carlo simulations were 
used to derive error estimates. 

(-H ' Results. Three transit light curves exhibit signatures of the occultations of dark spots by the planet during its passage across the 

stellar disk. The influence of stellar activity on transit depth is taken into account while determining system parameters. The radius of 
WASP-10 b is found to be no greater than 1.03+qoj Jupiter radii, a value significantly smaller than most previous studies indicate. We 
calculate interior structure models of the planet, assuming a two-layer structure with one homogeneous envelope atop a rock core. The 
' high value of the WASP-10 b's mean density allows one to consider the planet's internal structure including 270 to 450 Earth masses 

of heavy elements. Our new mid-transit times confirm that transit timing cannot be explained by a constant period if all literature data 
points are considered. They are consistent with the ephemeris assuming a periodic variation of transit timing. We show that possible 
starspot features affecting the transit's ingress or egress cannot reproduce variations in transit timing at the observed amplitude. 

Key words, stars: starspots - stars: individual: WASP-10 - planetary systems - planets and satellites: individual: WASP-10 b 
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^ ; 1 . Introduction and the mass of M„ = 0.75+°j* M Q (Johnson et al. [2009b . A ro- 

■»-4- . tational period of 1 1 .91 + 0.05 d was determined from photom- 
^T. , Exoplanets, whose orbits are aligned in such a way that they affected fe starspot . induced brightness modulation (Smith 

transit their host stars, constitute a group of worlds whose study gt a , 20Qg) The sim0ar iodicit of x L84 d was detected in 

O . can in turn help us to understand the exoplanetary internal radial velocities b Macie jewski et al. dMU, who re-analysed 

-H . structure and atmospheric properties (see e.g. Charbonneau et available trosc ic data These authors also found that the 

al. pEJ. The transiting exoplanet WASP-10 b (Christian et mhm eccentricit of the lanet is statistically indistinguishable 

al. 2009) was discovered by the SuperWASP survey (Pollacco from ^ ^ ^ ^ accQunt ^ Qt effect The 

et al. [20061). Us mass was initially measured to be 2.96*$ tem was found t0 be relatively young with an age of 27 ± 80 

Mj (Christian et al. |2009|) and then refined to 3.1 5ffjf Mj by Myr determined with the gyrochronology method (Maciejewski 



X 
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. Johnson et al. (20091 120101 ). Christian et al. ( 120091 ) reported the et a j 201 1 1 
planetary radius of 1.28^[!2g Rj and note that it is larger than 

interior models of irradiated giant planets predict. On the basis High-precision transit light curves of a planet moving across 

of a high-precision light curve obtained with the University of a s P otted stellar dlsk offer an outstanding opportunity to map 

Hawaii 2.2-m telescope, Johnson et al. (120091 derived a notice- the stars P ot distribution (Schneider |2000j Silva |2003J). When a 

ably smaller radius of 1.08 ± 0.02 flj. Further investigations by P lanet occults a starspot area on a star's surface, the flux in- 

Krejcova et al. (120101) and Dittmann et al. (t20T0b found a plane- creases and a characteristic feature in a light curve - a bump - 

tary radius close to the value reported by Christian et al. (120091 . 18 observed, as for example forHD 1 89733 b (Pont et al. |^TJ, 

The host star is a faint ( V = 12.7 mag) K5 dwarf located HD 209458 b (Sdva-Vaho 2008), GJJ214 b (Berta et al. [20111 

90 ± 20 pc from the Sun in the constellation of Pegasus. It has Carteret aLlMIl Kundurthy et al. 201 1 ), TrES-1 (Dittmann et 

the effective temperature of 4675 ± 100 K (Christian et al. 2009J al.[2009jl, CoRoT-2 b (Wolter et al. 2009; Huber et al.|20T(j), or 

WASP-4 b (Sanchis-Ojeda et al. 1201 II ). Unocculted dark spots 

* Based on observations collected at the Centro Astronomico outside the planetary path across the stellar disk may still be 

Hispano Aleman (CAHA), operated jointly by the Max-Planck Institut present on the stars's surface and affect the transit light curve, 

fur Astronomie and the Instituto de Astrofisica de Andalucia (CSIC). Compared to a spot-free case, spots reduce the effective stellar- 
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disk area, so a stellar radius is underestimated. The observed 
transits are deeper than a spot-free scenario predicts. This, in 
turn, results in the overestimate of the planet-to-star radius ratio 
(Carter et al. 1201 11 1. The opposite effect may be observed if facu- 
lae are considered. These active regions should introduce rather 
random variations not only in the transit depth, but also in the 
transit timing and duration. 

Interestingly, Maciejewski at al. (1201 II ) note that transit tim- 
ing of WASP-10 b cannot be explained by a constant period but 
by a periodic variation. The existence of a second planet that 
perturbs the orbital motion of WASP-10 b was postulated as an 
explanation of the observed transit-time-variation (TTV) signal. 
A solution with a perturbing planet of mass ~0. 1 Mj and orbital 
period of ~5.23 d, close to the 5:3 period commensurability, was 
found to be the most likely configuration. 

In this work, we present the analysis of four new high-quality 
light curves of WASP-10 b's transits. The aim of this work is to 
redetermine the system's parameters and verify the existence of 
the observed TTV signal. We also report on discovering starspot 
features in transit light curves and show the results of modelling 
the internal structure of the WASP-10 b planet. 

2. Observations and data reduction 

The 7?-band photometric monitoring of WASP-10 b's transits 
was performed with the 2.2-m telescope at the Calar Alto 
Observatory (Spain) during four observing runs in 2010 on 
August 3 and 6 and September 6 and 9 (programme H10-2.2- 
011). An additional four hours of monitoring during the out-of- 
transit phase were run on 2010 September 10 to check a photo- 
metric stability. During another two nights, which were granted 
to the project in 2010 November, poor weather conditions did 
not allow us to gather scientific data. The Calar Alto Faint Object 
Spectrograph (CAFOS) in imaging mode was used as a detector. 
It was equipped with the SITe CCD matrix (2048 x 2048, 24jt/m 
pixel, scale of 0.53 arcsec per pixel). The telescope was signif- 
icantly defocused and stellar profiles exhibited a doughnut-like 
shape. About 5 x 10 4 counts per second were recorded for the 
target, which allowed us to collect up to 3 x 10 6 counts in a 
single exposure. Applying defocusing was expected to minimise 
random and flat-fielding errors (e.g. Southworth et al. 12009b . A 
subframe limiting the effective field of view to the target and 
a nearby comparison star GSC 02752-00151 and 2x2 binning 
mode were used to shorten a readout overhead time to ~25 s. 
Including overheads, the observations were carried out with a 
cadence of between 66 and 83 s. Exposure times in individ- 
ual runs were kept fixed to avoid affecting transit timing. The 
colour index and brightness of the comparison star, located 3. '6 
to the east of WASP-10, were not far from those of the host 
star that minimised the photometric trends caused by the dif- 
ferential atmospheric extinction. The photometric stability of 
GSC 02752-00151 was verified by analysing its light curve pro- 
vided by the SuperWASP survey. We also checked its bright- 
ness against fainter stars available in the field and detected no 
variability. The stellar images were kept exactly at the same po- 
sition in the CCD matrix during each run thanks to auto guid- 
ing. Precise timing was assured by synchronising the computer's 
clock to Coordinated Universal Time (UTC) using Network 
Time Protocol software, which is more accurate than 0.1 s. The 
details of observations are presented in Table Q] 

Weather conditions during run 1 on 2010 August 3 were 
mainly photometric with occasionally thin clouds in the first half 
of the observing window. The second half of the run was affected 
by the light of the rising Moon. A portion of data was lost at 



the end of an egress phase because of technical problems with 
dome-slit segments. On 2010 August 6 (run 2), observations 
were performed in photometric conditions in dark time. The run 
was again interrupted by the same technical problems during the 
end of an ingress and the beginning of a flat-bottom phase. The 
subsequent two runs on 2010 September 6 and 9 were executed 
in photometric conditions with no moonlight contamination. On 
September 6, some exposures were lost at the end of the transit 
owing to a strip of clouds passing across the field. Exposures that 
were found to be affected by clouds were rejected. 

Our CCD frames were processed using a standard procedure 
including debiasing but not flat-fielding, which does not improve 
measurements in the case of defocused and auto-guided observa- 
tions and may also degrade the data quality (see e.g. Southworth 
et al. 120101) . Our tests indeed showed that dividing by a flat-field 
frame had no detectable effect on a final light curve. The stellar 
flux was spread over a ring of at least 20 binned pixels (21'.'2) 
in a diameter. The magnitudes of WASP-10 and the compari- 
son star were determined with differential aperture photometry. 
After analysing a wide range of apertures, we found that begin- 
ning from the aperture radius of 9 binned pixels, the photometric 
scatter exhibited a flat minimum (or rather a plateau), which is 
in practice insensitive to greater aperture sizes. The per-exposure 
photometric uncertainties in the individual measurements were 
in a range between 0.59 and 0.82 mmag. A visual inspection 
of the Digitized Sky Survey (DSS) images revealed at least two 
faint neighbour stars located 12'.'5 and 1579 away from WASP- 
10. Their photographic F-band magnitudes are 19.6 and 19.9 
mag (Zacharias et al. 120041) . hence negligible relative to those of 
WASP-10, whose brightness in the same band is 11.8 mag. To 
minimise the neighbour-star contamination, the light curves with 
minimal apertures for which the lowest scatter was achieved, i.e. 
9 binned pixels (1076), were taken as the final ones. The light 
curves exhibited slight baseline gradients caused mainly by the 
differential atmospheric extinction. They were approximated by 
the first or second order polynomials fitted to out-of-transit mea- 
surements and then subtracted from individual light curves. This 
method also removed a possible brightness gradient caused by 
the host-star variability and the possible effects of moonlight. 
However, taking into account the small amplitude of the stellar 
variability (up to ~20 mmag, Maciejewski et al. 1201 11 see also 
Section 3.2) and its relatively long timescale of ~12 d (Smith 
et al. 120091 ), its contribution remains negligible on the timescale 
of the observing windows. To estimate the quality of the light 
curves, the root-mean-square (rms) of data points at the out-of- 
transit phase was calculated for individual light curves (Table[TJ. 

3. Results 

3. 1 . Light-curve modelling 

A visual inspection of the first three transit light curves (Fig.Q]) 
reveals features that affect the shape of transits and could be at- 
tributable to starspot features. The data from the fourth run were 
found to be of the highest quality exhibiting the lowest data 
scatter at the out-of-transit phase. This was achieved thanks to 
photometric conditions during the run and long exposure time. 
Moreover, these data were found to be almost free of any de- 
formations, thus were used to perform the initial modelling. 
We used the JKTEBOP code (Southworth et al. 2004a 2004b), 
which is based on the EBOP programme (Eclipsing Binary Orbit 
Program; Etzel ll98"Tl Popper & Etzel 1981). The initial param- 
eters of the system were taken from Johnson et al. ([2009). The 
best-fit model was found by applying the Levenberg-Marquardt 
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Table 1. The summary of observing runs. 



Run 


Date UT 


Time UT 


exp 


X 


^exp (SJ 


Trad (S) 


FoV 


rms 00 t (mmag) 


1 


2010 August 03 


21:44-01:48 


165 (175) 


1.94 -» 1.03 


40 


67 


7:83x4^88 


1.32 


2 


2010 August 06 


23:28-03:40 


213 (217) 


1.24 -> 1.01 -> 1.04 


40 


66 


7:83 x 4:88 


1.16 


3 


2010 September 06 


22:20 - 02:04 


164 (189) 


1.11 -> 1.01 -> 1.08 


50 


71 


5:57 x 334 


1.06 


4 


2010 September 09 


00:12-04:26 


184 (185) 


1.01 -> 1.64 


60 


83 


6:02 x 3:73 


0.79 



Notes. N cxp - the number of useful scientific exposures, a total number of exposures in parentheses, X - airmass changes during the run, r exp 
- exposure times, T cad - the cadence, FoV - the field of view, dates in UT at the beginning of the night, rms oot - the data scatter during the 
out-of-transit phase per point. 




-0.08 -0.04 0.00 0.04 0.08 -0.08 -0.04 0.00 0.04 0.08 

Time from mid-transit (d) 

Fig. 1. Light curves for four transits of WASP-10 b with models based on the solution obtained for run 4 data and employing the 
quadratic limb darkening law. The residuals are shown in bottom plots. 



least squares procedure. To derive the parameter errors, we used 
procedures implemented in the JKTEBOP code and described 
in detail by Southworth (120081) . Firstly, we ran 10000 Monte 
Carlo (MC) simulations to account for Poisson noise with a 
bootstrap approach. Secondly, we employed the residual-shift 
(prayer-bead) method (e.g. Desert et al. 1201 II Winn et al.[2009) 
to check whether a light curve is affected by correlated noise. A 
spread range in a given parameter of within 68.3% was taken as 
its lcr error estimate. 

Six parameters defining the transit shape were fitted directly: 
the sum of fractional radii of the host star and planet r, + r^, 
the ratio of these quantities k = r\,/r*, the orbital inclination i, 
the mid-transit time To, and the linear u and non-linear v limb- 
darkening coefficients (LDCs). The fractional radii are defined 
as r„ — — and n, = — , where R t and Rv. are the absolute 
radii of the star and the planet, respectively, and a is the orbital 
se mi-major ax is. A value a = 0.03781*o^^ AU (Johnson et 
al. 2009 2010) was used in subsequent calculations. Combining 
i and R* allowed calculation of the transit parameter b — cos i. 
The mid-transit times were transformed from Julian Date (JD) 
based on UTC into Barycentric Julian Date (BJD) based on 



Barycentric Dynamical Time (TDB) using the on-line convertetQ 
by Eastman et al. ( 120101 . 

We checked linear, quadratic, logarithmic, and square-root 
limb-darkening (LD) laws and various configurations of LDCs, 
for which theoretical values were bilinearly interpolated from 
tables by Van Hamme ( 1993 ). We considered LDCs kept fixed 
during fitting, u as a free parameter while v kept fixed, and al- 
lowed both LDCs to vary. The best-fit model could be obtained 
using the square-root LD law with both LDCs being fitted, but 
the result turned out to be unphysical (total limb darkening at 
the limb of the star was greater than 1). Keeping v fixed at a 
theoretical value for the square-root LD law resulted in a poorer 
goodness of the fit than applying the quadratic LD law for which 
both LDCs were allowed to vary. Therefore, we used the lat- 
ter configuration to reproduce the light curve. The MC errors 
were found to be up to 50% smaller than those estimated by 
the prayer-bead method. This indicates that the light curve is af- 
fected not only by Poisson noise but also by additional correlated 
noise. Therefore, our prayer-bead error estimates were taken as 
our final errors. The sub-millimagnitude precision was achieved 



1 http://astroutils.astronomy.ohio-state.edu/time/utc2bjd.html 
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Fig. 2. A light curve acquired for WASP-10 during an additional 
run without transit. The data were processed in a similar way as 
for transit light curves. The 55-s exposure times were used. 



with the rms scatter of 0.74 mmag per point. We note that de- 
rived LDCs have reasonable values that are consistent to within 
lcr with theoretical values interpolated e.g. from Claret (2000) 
tables for /?c-band filter and a K-type dwarf. Fitting LDCs al- 
lows us to account for their contribution to the error budget. 
Southworth (2008) shows that keeping LDCs at theoretical val- 
ues can underestimate the errors of other parameters. 

The best-fit model based on data from the fourth observing 
run was used as a template and fitted to remaining light curves 
for which only the mid-transit time was allowed to vary. The 
residuals are plotted in Fig.Q] Both transits observed in August 
exhibit flux bumps that may be the results of occulting dark- 
spot areas on the stellar surface by the planetary disk. The neg- 
ative residuals in run 2 data suggest that in this case the transit's 
depth is underestimated by the run 4 model. The opposite sit- 
uation is observed in the light curve acquired on September 6 
where, besides the starspot signature, the transit depth tends to 
be smaller than the run 4 model predicts. Figure|2]shows the light 
curve obtained for WASP-10 on September 10 during an out-of- 
transit phase. During this almost four-hour-long run, the rms of 
0.81 mmag was achieved. The residuals show no features that 
could be attributed to starspots. As this run was also affected by 
clouds, the relatively far more satisfactory photometric stability 
strengthens the reliability of starspot detection in runs 1-3. We 
also exclude observing conditions or technical problems as the 
cause of the observed bumps. 

Variable transit depth may be caused by apparent changes 
in the effective stellar surface (hence fitted fractional stellar ra- 
dius), which varies because of the dark-spot coverage. To inves- 
tigate this scenario, we cut portions of the light curves exhibiting 
spot signatures and then repeatedly fitted models and estimated 
the errors. Constraints on the limb darkening were applied be- 
cause the quality of the reanalysed three light curves is slightly 
lower than for run 4 data and data sets became significantly in- 
complete after cutting. We considered only the quadratic LD law 
and LDCs were kept fixed at values from the run 4 model. Their 
contribution to the error budget was taken into account by per- 
turbing them by ±0.1 around fixed values with the assumption 
of a flat distribution. Figure [3] shows individual best-fit models 
with residuals, which, if uncut light curves are considered, con- 
tain the distilled spot features. Results of the fitting procedure 
run for individual light curves are collected in Table [2] 

3.2. Physical properties of the system 

The highest-quality light curve of run 4 allowed us to find the 
best-fit LD law and its LDCs but other system parameters (in 
particular, the planet-to-star radius ratio k) may be affected by 
stellar activity. Although no spot signatures are visible during 
this transit, spot complexes may exist on the stellar surface out- 
side the projected path of WASP-10 b. To check this scenario, the 
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Fig. 3. Light curves for three transits of WASP-10 b with individ- 
ually fitted model light curves after rejecting data points affected 
by stellar spots (grey points). The residuals that were used to 
model starspot features are plotted in the middle graphs. In the 
bottom graphs, the final residuals are plotted. 



average brightness of WASP-10 outside transits was determined 
in individual nights covered by our observations. Figure|4]shows 
the light curve that clearly exhibits data-point scatter with the 
peak-to-peak amplitude of ~20 mmag. These brightness varia- 
tions can be closely approximated with a sinusoid whose pe- 
riod is equal to the star rotation period and the peak-to-peak 
amplitude is 28 + 13 mmag. Individual magnitude errors (stan- 
dard deviations) were taken as weights while running the fitting 
procedure. The derived amplitude is comparable to the greatest 
value reported so far, i.e. to 20.8 +1.3 mmag observed for the 
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Table 2. Parameters of transit light-curve modelling. 



Parameter 


run 1 (Aug. 03) 


run z (Aug. Uo) 


run j (oep. Uo) 


run 4 (oep. UV) 


Sum ol stellar and planetary fractional radii, r, + r\ 


n nntt+ll 0030 

1 SQfi +U * UU3U 

u,1J7U -0.0029 


°- 0965 -fissf 




0973+° mli 

Q QQ13 
n 1 ^OA+0 00 '8 
U.13SO_ 00027 


Planet-to-star radius ratio, k 


, ,^ + (1 (1(172. 

u ' iUU/ -0.0032 


1 ShV umu 
u,UUJ -0.0022 


Stellar fractional radius, r. 


0.0833!™** 




0-0826™ 


0.0840+°;°°;? 


Planetary fractional radius, rt 


n n 1 ti +0 0007 
u,uu -'-0.0003 


o ni ^4 


n ni?o+ 00006 

\).\)iz.y -0.0001 


j-i n 1 n i +0 0004 
01 T^ +UUUU4 


Orbital inclination, i (deg) 


89.0 " 


QQ T+O.? 


89.5 " 


88.9 „, 


Transit parameter, b (R,) 


0.19*"« 


1 ^+0 14 

() 1 n ' 


. n i OS 
1 ()+ UUB 


n <^*-»+0 0Q 

-0.08 


Transit duration" , d lT (min) 


115.9+JJ 


1 16.3+42 


116.5*|| 


ii6.r„ 


Ingress/egress duration* , d w (min) 


19.2+H 


1 9 2+ 3 - 5 


18.4+2-3 


19.4+| 3 L 

n £c+0~07 
°- 65 -0O9 


Linear LDC, u 


0.65 r 


0.65 r 


0.65 c 


Non-linear LDC, v 


0.16 1 ' 


0.1 6 C 


0.1 6 C 


16+°- 6 

49 S85Q 8 +0 0008 
— Q.QQQQ9 


Mid-transit time, To (JDutc) 


12.47488+°;°™° 


15 56763+ - 00015 


4(S 4Q^47+0-00007 
46.49882« 


Mid-transit time, T (BJD TD b) 


12.47861+™}° 


15 57156+™' 4 

' iJU -0.00013 


49 59 140+0 00088 


Epoch rf 


341 


342 


352 


353 


Timing residual, O-C (d) 


-0.00009 


+0.00013 


-0.00009 


-0.00006 


rms (mmag) 


1.42 


1.11 


0.97 


0.74 


X red 


3.40 


1.99 


2.06 


1.53 



Notes. The mid-transit times, T , are given as 2455400+ . No spot feature cutting was applied for run 4. 

<n> Time between the middle of an ingress and the middle of an egress, defined as d u = r,P^ 1 Vl - b 2 , where is the orbital period of 
WASP-10 b taken from Maciejewski at al. d201 It . (i) Time between the first and second or third and fourth contacts for ingress or egress, 
respectively, defined as d a = ^-^t . where P\, is the planet orbital period. <r) Value adopted from the run 4 fit, permuted by ±0.1 on a flat 
distribution. (rf) Numbering according to the ephemeris given by Christian et al. (2009 I. 



Epoch 



340 345 350 355 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 j 




■ I ■■■■ I ■■■■ I .... I I . ■ 

10 20 30 40 50 

BJD TnR -2455400 



Fig. 4. The average brightness of WASP-10 at an out-of-transit 
phase relative to the comparison star GSC 02752-00151. The 
continuous line sketches the sinusoidal variation caused by the 
rotational modulation due to spots. Individual runs are marked 
with 1-4 numbers, while the additional one during the out-of- 
transit phase is labelled with A. 



2004 observing season (Maciejewski at al. 1201 11 1. Run 4 was 
observed during a phase at which WASP-10's brightness was re- 
duced as a result of new starspots appearing on the visible stellar 
hemisphere. This finding indicates that k derived from these data 
must be overestimated. The run 1 and 3 data show that even dur- 
ing bright phases, spot features are visible in transit light curves. 
This finding suggests that the purely photospheric surface was 
observed in neither run. The photospheric brightness level can- 
not be determined, so the reliable rescaling of light curves cannot 
be performed (Czesla et al.[2009 ). We note that the LDCs may be 
affected by stellar activity, e.g. faculae that produce limb bright- 
ening and, if numerous, may have a non-negligible contribution 
to the overall limb darkening. 



The presence of starspots located outside the transit chord 
leads to an overestimate of the transit depth and hence the planet- 
to-star radius ratio (e.g. Czesla et al. 2009). Although system 
parameters are consistent to within l<x between individual runs, 
the smallest value of k was obtained for run 3. The star was then 
at a bright phase (Fig. 0), so a minimal fraction of its surface 
must have been covered by spots. The run 3 results were chosen 
as the most representative of the WASP-10 system. They were 
also expected to be the closest to the real values. Therefore, we 
used them to redetermine the planetary, stellar, and geometrical 
properties. Our final results are collected in Table [3] 

Our parameters of the WASP-10 system differ by 0.6-1.9cr 
from those reported by Johnson et al. (2009 2010) with only 
the planet gravitational acceleration value differing by more than 
2<x. The only exception is whose value by Johnson et al. 

(20091 120101 ) is ~ V2 times greater than our determination. The 
reason for this discrepancy is the factor of V2 missing in the 
formula defining the equilibrium temperature in Johnson et al. 
(2009). The high-quality light curve analysed by Johnson et al. 
(2009) exhibits no signs of spot occultations and, as an extrap- 
olation of the brightness modulation in Fig. [4] suggests, was ac- 
quired during a bright phase, i.e. when the spot coverage was 
at its smallest. System parameters reported by Christian et al. 
([20091) and Krejcova et al. d20T0| ) differ by 1.3-5.4cr from those 
derived by us. The planet-to-star radius ratio and star radius have 
the greatest impact on other parameters. The light curve anal- 
ysed by Dittmann et al. d20T0T >, whose k = 0.16754 + 0.0 0060 
was found to be close to values of Christian et al. (2009J and 
Krejcova et al. (120101 ). exhibits some features attributable to 
starspots. Studies by Christian et al. (2009) and Krejcova et 
al. ( 120101 ) based on average light curves that blurred the influ- 
ence of activity-induced features. Although our radius estimate 
of Rb = 1.03+QQ3 Rj is the smallest determined so far, it is de 
facto the maximal limit because we could take into account only 
a fraction of the whole spot contribution. If the greatest value of 
k were considered (i.e. for the run 2 data set), the planetary radius 
would be greater by 3%, i.e. well within error bars. Hence, the 
difference seems to be negligible in practise. Using the Roche 
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Table 3. Physical properties of the WASP-10 system derived from light-curve modelling. Literature determinations are cited for 
comparison. 



Parameter 




This work 


Chustian et al. (zU09) 


Jonnson et al. (zUUV, 2U1U) 


Krejcova et al. (20101 


Planet radius, Rb (/Cj) 




1 m+0 07 


i 90+0.077 


1.08 ± 0.02 


1.22 ± 0.05 


Planetary mean density, pb (pj) 




r\ a +0 4ft 

2 -ns 

1 QQ+0.04 

■J.OB „ fl , 


1 - 4 -'-0.29 


2.35 ± 0.15 




Planet gravitational acceleration, log 


gb (cgs) 


J. 02 ± U.Uo 


o.ozo ± U.Ulz 




Equilibrium temperature, T' (K) 






1119^|| 


1370 ± 50 




Orbital inclination, i (deg) 




89.5 + °4 


86.9 +06 


88.49 + ^ 2 


87.3 ±0.1 


Transit parameter, b (R,) 




10 +oi)8 


0.568+ ; jj 4 


n 9QQ+0.020 




Planet-to-star radius ratio, Rt/R, 




1565+ 00033 


0.170 ±0.002 


1 59 18 +0 00050 


0.168 ±0.001 


Star radius, R, (R B ) 




67 +0 03 


n 77C+0.043 

-0.040 


0.698 ± 0.012 


0.75 ± 0.03 


Scaled semi-major axis, a/R, 




12 H+° 24 


10 23 +(W0 

"•"-492 


1 1 65 +009 


10.64 ±0.12 


Mean star density, p, (p Q ) 




2 48 +a26 


1 51 +025 

-Q-2Q 


2.20 ± 0.063 




Star gravitational acceleration, log g, 


. (cgs) 


4 66 +0 06 


4 < 1 +0.06 


4 fi 97+o.oioi 





Notes. The equilibrium temperature, 7"' , was derived by assuming the effective temperature of the host star Tgg = 4675 ± 100 K (Christian et al. 
20091 and simplified relation T = r cff V''*/2 (Southworth 20101. Some quantities that were needed for calculations and could not be determined 
from the light curve analysis (e.g. planetary mass, semi-major axis) were taken from Johnson et al. d2009l[20T0t . 



model, Budaj d20TTb has found that the shape of WASP-10 b is 
not far from the sphere. The ratio of the planetary radii at the 
sub-stellar point to that at the rotation pole is 1.00132. Thus, the 
planet radius determined from the transits is practically the same 
as the effective radius which is defined as the radius of the sphere 
with the same volume as the Roche surface used in theoretical 
calculations. Hartman (2010) discovered a correlation between 
the surface gravitational acceleration of hot Jupiters and the ac- 
tivity of host stars. According to this relation, the surface gravity 
of WASP-10 b, loggb = ^ (cgs), predicts a strong chro- 

mospheric activity of the WASP-10 star with logR' HK > -4.3. 
The star is expected to have pronounced emission features in the 
H and K Ca II lines that make it an interesting target for spectro- 
scopic studies. 

The satisfactory consistency is obtained when the derived 
gravitational acceleration of the host star, logg, — 4.64^^, is 
compared to theoretical predictions. We interpolated the solar- 
metallicity isochrones by Girardi et al. (2002 ) to the age of the 
WASP-10 system and derived a value of ~4.69, which is within 
lcr error bars of the value reported in this work. The value de- 
rived by Johnson et al. (20091 120101 ) is close to the theoretical 
one, while the result brought by Christian et al. (2009 ) is signif- 
icantly underestimated. 

3.3. Internal structure 

We calculate interior structure models of WASP-10 b assuming a 
two-layer structure with one homogeneous envelope atop a rock 
core. We aim to determine the resulting possible mass of heavy 
elements in both the core and the envelope as a function of the 
planet radius, which we vary between the small value of 1 .03 Rj 
(this work) and 1.28 Rj found by Christian et al. (120091), while 
the mass is kept invariant at 2.96 Mj (Christian et al. 2009). 

For H and He we use the interpolated equation of state (EOS) 
by Saumon et al. (1995 ), for heavy elements in the envelope we 
scale the He EOS by a factor of four in density, and for rocks 
in the core we use the rock EOS by Hubbard & Marley ( 1989). 
This is a choice of convenience, as we are interested in the re- 
sponse of the planet's metallicity to the planet's radius more than 
to achieve a high accuracy in these numbers for a single model. 
The composition of heavy elements is essentially unknown, and 
a different choice of the EOS of metals can affect the calculated 
mass of metals by ~ 30%. 



An interior model of a close-in exoplanet cannot stand alone 
without a model for the irradiated model atmosphere, as the at- 
mosphere determines the entropy of a gas giant's interior. The 
hotter the atmosphere, the larger the entropy, and the more 
heavy elements are required to meet the observed transit radius. 
Common features of the atmospheric pressure-temperature pro- 
file of evolved hot Jupiters around Sun-like stars are the presence 
of a temperature inversion at high altitudes at sub-mbar pres- 
sures, the sub-adiabatic temperature profile in the radiative re- 
gion where the transit radius is observed (Fortney et al. 2003 ), 
the development of an isothermal region around 1-100 bars de- 
pending on an age and orbital distance, and finally the onset 
of the adiabatic, convective interior, which extends down to the 
core in the likely absence of internal layer boundaries in hot exo- 
planets. We use Fig. 3 in Fortney et al. (2007 ) to construct an at- 
mosphere profile by interpolation, where we increase the orbital 
distance from a = 0.037 AU to 0.07 AU in order to conserve 
the total energy flux received by WASP-10 b from its parent K5 
star. The resulting isothermal region of our atmosphere profile 
is at 1400 K and begins at 2 bars. As the isothermal region ex- 
tends deeper into the planet with increasing planet age (Fortney 
et al. 120071 ). we expect a shallow or even absent isothermal re- 
gion for the young WASP-10 b. We use Fig. 2 in Fortney et al. 
(120071 ) to estimate the onset of the adiabatic interior, which we 
parametrise by a pressure level P at j. The equilibrium temperature 
of WASP-10 b is 850-980 K for albedos between 0.3 and zero. 
Hence, an effective temperature of 1000 K is certainly a lower 
limit for WASP-10 b. For our choice of planet radii, the surface 
gravity is loggh = 3.62 (i.e. gb — 42 m s~ 2 from Christian et 
al. 12009) and loggb — 3.88 (i.e. gb, — 76 m s~ 2 determined by 
us). For a young hot Jupiter with logg = 3.6 (g — 40 m s~ 2 ) and 
r eff = 1000 K at 0.1 AU, Fortney et al. ( 120071 predict P ad * 3 
bar. We vary P a + from 0.5 to 10 bars to cover the uncertainty re- 
lated to a possible higher and larger gravity. The procedure 
for calculating the internal profile and the core mass is the same 
as described in Nettelmann et al. d201 II ). The envelope heavy 
element mass fraction is set to the solar value of 1.5%, which 
also agrees with the stellar metallicity to within the error bars. 
This choice gives the maximum core mass as a higher envelope 
metallicity reduces the core mass when all other parameters are 
not varied. 

Figure|5]shows the results of our calculations. The total mass 
of metals Mz of 270-450 M ffi is needed to reproduce the radius 
value reported in this work. This large amount is challenged by 
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Pad [bar] 



Fig. 5. Core mass (solid lines) and total mass of metals (dashed 
lines) of WASP- 10b interior models with solar metallicity enve- 
lope, both in M & . The models differ in the onset of the adiabatic 
interior as parametrised by the pressure level P a d- The observed 
transit radius of this work requires to increases the total mass 
of metals by ~300 M e compared to the larger radius from the 
discovery paper. 



the core-accretion formation where the maximum mass of avail- 
able heavy elements is 80 M @ for the parameters of WASP-10 b 
(Leconte et al. 2010). In contrast, the larger radius scenario pre- 
dicts amass of heavy elements that is smaller by ~300M e , which 
corresponds to a heavy element enrichment of (12 + 11) times a 
stellar value, i.e. more in line with Jupiter. 

3.4. Spot feature modelling 

The spot-occultation features (Fig. [3]) may be approximated with 
a simplified transit-like model assuming a flux bump instead of 
a flux drop. To model these features, an approach similar to that 
presented by Kundurthy et al. (201 1) was used. The bumps were 
modelled with the occultsmall routine of Mandel & Agol ( 2002) 
and the Levenberg-Marquardt non-linear least squares fitting al- 
gorithm provided by the Exoplanet Transit Databas^l (Poddany 
et al. 120101 . In this case, the limb darkening was assumed to be 
zero and the mid-occultation time, occultation height, and dura- 
tion counted from the begin to the end of the phenomenon were 
fitted. The error estimates of these parameters were taken from 
a covariance matrix. The best-fitting models are plotted in Fig. [3] 
(middle plots) and occultation parameters are printed in Table[4] 
The final residuals are shown in the bottom plots of Fig. [3] The 
height of the bumps was found to be between 2 and 3 mmag. The 
durations of the events are slightly shorter than a sum of tran- 
sit ingress and egress durations equal to ~40 min (see Table [2]). 
Moreover, no flat top phase is visible. These findings are most 
consistent with the grazing scenario and non-central occultations 
of single spots whose sizes are comparable to or larger than the 
planetary disk projected onto the stellar surface. The occulted 
starspots were found to be located near the centre of the stellar 
disk at longitudes within 20° of the line of sight. The height of 
the bumps was determined at the 5-9-cr level. The significance 
of these bumps can be estimated from the ratio of their height 
to the final light-curve scatter represented by the rms. The de- 
tections of starspot signatures in run 1, 2, and 3 data were found 
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at the 1.4-, 2.9-, and 2.0-cr level. The run 2 event was the most 
pronounced while the bump in the run 1 light curve was affected 
by lower quality data. While reducing data, we excluded instru- 
mental effects being responsible for the observed features in the 
light curves. 

The flux bumps can be translated into the area of a spot oc- 
culted by the planetary disk. Assuming that the spot contrast is 
equal to 0.7, a value typical of the Sun, the occulted spot area was 
found to be in a range between 0.6% and 0.9% of the stellar sur- 
face (or 25% and 37% of the planetary disk). These inhomogen- 
ities alone cannot produce the observed rotational modulation in 
photometry. To reproduce the observed amplitude of ~28 mmag, 
a distinct spot of a size of either ~9% or ~2.6% of the stellar 
disk would be needed for a solar contrast or completely dark 
spot, respectively. As the amplitude of the rotational modulation 
was found to be one of the greatest observed for WASP-10 so 
far, the radial velocity variations with a semi-amplitude much 
greater than ~70 m s _1 , a value reported by Maciejewski et al. 
(1201 11 1, would be expected if simultaneous spectroscopic obser- 
vations were done. The starspot features observed in the transit 
light curves are consistent with the observed radial velocity scat- 
ter and the photometric variability of the host star. 

3.5. Transit timing 

Four new mid-transit times in conjunction with the literature data 
allowed us to refine the linear ephemeris. We obtained the mid- 
transit time for initial epoch T = 24554357.85790 + 0.00010 
BJDtdb and the orbital period P h = 3.09272963 ± 0.00000035 
d. Individual timing errors were taken as weights while running 
the linear regression. Our four new mid-transit times could be 
used to verify the TTV signal reported by Maciejewski et al. 
( 1201 II) . Figure [6] shows the observation minus calculation (O- 
C) diagram plotted for all data according to the new ephemeris. 
Although our four new points are located close to the lin- 
ear ephemeris, they also agree with the ephemeris given by 
Maciejewski at al. (201 1 ). The remaining two runs, which were 
granted to us and then lost because of poor weather conditions, 
had been expected to cover epochs during opposite phases of the 
TTV signal and could have been used to determine precisely its 
amplitude. It is worth noting that all four transits happened ~5 
min later than the linear ephemeris by Dittmann et al. {2010) 
predicts. This finding confirms the TTV signal of WASP-10 b 
and is consistent with the model proposed by Maciejewski et al. 
( 120111) . 

To estimate the influence of starspot occultations on transit 
timing, the run 4 light curve was artificially deformed by inject- 
ing a flux bumps randomly located in the time of transit. The du- 
ration of artificial occultations was assumed to be ~40 min and 
its height 3.07 mmag, i.e. values of the most prominent feature 
observed in run 2. Two separate scenarios were considered: an 
affected first half of a flat bottom and ingress/egress. In the first 
case, 2000 deformed light curves were produced and then fitting 
procedure was performed. The sum of stellar and planetary frac- 
tional radii, planet-to-star radius ratio, orbital inclination, and 
mid-transit time were allowed to vary while the quadratic LDCs 
were kept fixed at values from the run 4 model. The mean dif- 
ference in To, 6Tq, was found to be negligible with a value of 
0.00006 + 0.00004 d (i.e. 0.09 ± 0.06 min), where the error is 
the sample standard deviation. This finding was also confirmed 
by our analysis of run 1-3 light curves: the mid-transit times 
derived from original and starspot-feature-cut light curves were 
found to be similar well within lcr error bars. The planet-to-star 
radius ratio was found to be significantly smaller by about 5cr. 
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Table 4. Spot-occultation parameters. 



Parameter run 1 (Aug. 03) run 2 (Aug. 06) run 3 (Sep. 06) 

Height (mmag) 1.99 ± 0.44 3.07 ± 0.32 1.87 ±0.31 

Duration (min) 34.3 ± 7.9 37.2 ±4.0 31.6 ±5.5 

Mid-occultation time - T (d) +0.0008 ± 0.0014 -0.0074 ± 0.0007 -0.0141 ± 0.0010 

Percentage of stellar disk area" 0.60 ±0.13 0.9 ±0.1 0.6 ±0.1 

Longitude* (°) +1.1 -10.5 -20.3 

rms (mmag) L43 L07 035 

(o) The spot contrast of 0.7 (typical for the Sun) was assumed. ib) Defined as -90° at ingress, 0° at mid-transit time, and +90° at egress. 



This result meets expectations because injected bumps increase 
the averaged brightness during the flat bottom phase, hence de- 
crease the averaged transit depth. 

The procedure was repeated for ingress phase based on 1000 
deformed light curves and yielded 6T = 0.00062 + 0.00005 d 
(i.e. 0.89 ± 0.07 min), indicating that the apparent mid-transit 
times happened noticeably later than the spotless ephemeris pre- 
dicted. The orbital inclination was found to be underestimated 
leading to a greater value of the transit parameter. Moreover, 
transits appeared to be shorter on average by 1.5 + 0.6 min. 
Because of the symmetry of the phenomenon, the same results 
can be reached for the deformed second part of the flat bottom 
and egress. In these cases, the apparent mid-transit times would 
happen earlier. This phenomenological approach gives the upper 
limits to any possible deviations in timing. As Sanchis-Ojeda et 
al. ( 1201 11 1 show, the influence of occulted starspots located close 
to the stellar limb is smaller by a factor of 3-5 because of the 
effects of limb darkening and geometrical foreshortening. 

4. Discussion 

Three of our four transit light curves exhibit signatures of 
starspots, while Maciejewski et al. d201 II) detected no variation 
in the transit shapes of nine light curves. Most of their data sets 
have a data scatter of around or larger than 2 mmag. Such a pre- 
cision seems to be insufficient for detecting flux bumps reported 
by us. Their only light curve with the data scatter comparable to 
ours exhibits no features that could be related to an occultation 
of an active region. However, this finding should be treated with 
caution as a final part of the egress and out-of-transit phase was 
not covered by observations. 

In the case of defocused images, the light contamination 
from nearby faint stars may make a transit shallower and hence 
lower the planet-to-star radius ratio. In the neighbourhood of 
WASP-10, there are two faint stars that could contribute light to 
the WASP-10's flux (see Section 2). Assuming that the total flux 
of these neighbouring objects (based on their photographic mag- 
nitudes) contributed to the WASP-10 light curve, the planetary 
radius would be underestimated by 0.07% - a value negligible 
in practice. 

The high mean density of WASP-10 b (2.94+° f 5 p 3 ) in con- 
junction with the relatively young age of the WASP-10 sys- 
tem testifies for the planet's heavy-element core of a signifi- 
cant mass, as discussed by Johnson et al. (120091) . The mod- 
els of internal structure by Fortney et al. (2007) and Baraffe 
et al. ( 120081) predict a substantial core of the mass >100 M ffi 
for WASP-10 b. In the mass-radius plot (Fig. 0, WASP-10 b 
is located in a sparsely populated area covered by deflated exo- 
planets. Our models of the WASP-10 b's internal structure show 
that for the small radius reported in this work, the total mass 
contained in metals (270^450 M ffi ) is much higher than that for 
Jupiter (12-A1 M ffl ; Saumon & Guillot 2004) and in terms of 
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Fig. 7. The mass-radius diagram for currently known transiting 
exoplanets (grey dots). The black spot marks WASP-10 b. The 
density contours of Jupiter (upper) and Earth (lower) are plotted 
with blue lines. The red line sketches the planetary isochrone 
for a 100 M ffi core planet, the interpolated equivalent semi-major 
axis of WASP- 1 b, and the age of 300 Myr (Fortney et al. l20U7l l. 



planet mass fraction (M z /M p i anet = 0.28 - 0.47), WASP-10 b 
is located between Jupite r (0.04-0.12) and HD149026 b (0.43- 
0.72, Ikoma et al. 2006). Since Mz is quite insensitive to the 
assumed distribution of heavy elements between the core and 
envelope, the Mz values can be used to derive an overall en- 
richment of 35 + 25 times the stellar metallicity Z* for which 
[M/H]=0.02 + 0.3 (Christian et al. 120091 . The extremely high 
core mass of > 250 M e indicates that the envelope is not of solar 
metallicity but of a few times Z Q , at least. 

If the spins of both the stellar rotation and planetary orbital 
motion are roughly parallel to each other, this configuration may 
allow one to either determine or refine a period of stellar rota- 
tion by analysing the starspot distribution probed by a transiting 
pla net (S ilva-Valio 120081 Silva-Valio & Lanza I20TT1 Huber et 
al. 1201 0). In a significant fraction of extrasolar systems, both 
axes were found to be misaligned (see e.g. Hebrard et al. 1201 11 ) 
and this could also apply to the WASP-10 system. The ratio of 
periods of WASP- 10 b's orbital motion to stellar rotation is close 
to 1:3.8, hence the occultations of the same spot near the centre 
of the stellar disk (as in runs 1-3) are impossible to observe dur- 
ing two consecutive transits. However, assuming that the active 
regions exist longer than a few stellar rotation periods and the 
system is not significantly tilted, we found that features observed 
during runs 1 and 3 (the time span of about three rotation peri- 
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Fig. 6. The O-C diagram for mid-transit times of WASP-10 b. The mid-transits published in the literature are marked with open 
symbols. The filled ones denotes new mid-transit times determined by us. The periodic signal predicted by Maciejewski at al. (1201 II) 
is sketched with a continuous line. The dashed line shows the refined linear ephemeris and the dotted one marks the linear ephemeris 
provided by Dittmann et al. ( 120 101) . 



ods) could be attributable to the same spot complex. In this case, 
mid-occultation times would refine the rotation period of WASP- 
10 to 1 1.85 d and 10.87 d for prograde and retrograde orbits, re- 
spectively. The former value agrees with the period known from 
photometric and spectroscopic observations. 

About 20% of the stellar disk is probed during transits. If the 
system is aligned, the transit chord covers a latitudinal band at 
the stellar equator (up to about +10° in latitude). Systematic and 
long-timescale studies of the transit profile morphology would 
allow one to investigate a spot distribution during the cycle of 
stellar activity. 

A typical full width at half-maximum for WASP-10 cross- 
correlation functions (CCFs) obtained with the SOFIE spec- 
trograph in its high-efficiency (HE) observing mode by cross- 
correlating the WASP-10 spectra with a K5V template was found 
to be in a range of 10.2-10.4 km s" 1 (Christian et al. 2009). 
This translates into the mean projected stellar rotational veloc- 
ity vsin;» = 4.1 ± 1.0 km s _1 . Here we used the calibration of 
the projected rotational velocity for the HE observing mode of 
the SOPHIE spectrograph by Boisse et al. (2010) with a typical 
error of 1.0 km s . We also assumed the star's intrinsic colour 
index to be (B - V)o = 1.15 + 0.10 mag, a value derived from 
the spectral type and Schmidt-Kaler (1982) calibration for main 
sequence stars. Here the formal error estimate reflects the un- 
certainty in one spectral class. We also checked the consistency 
of the WASP-10's spectral type with its near-infrared colour in- 
di ces fro m the 2-Micron All Sky Survey (2MASS, Skrutskie et 
al. 12006 ) using the relation given by Ducati et al. (1200 II ). Using 
the rotation period determined by Smith et al. (2009), the mini- 
mal radius of the star appears to be 0.98^Q 3g R Q . This value is by 
40% greater than our photometric analysis predicts. Both val- 
ues are marginally consistent within the error bars. Additional 
spectroscopic observations are needed to remove this discrep- 
ancy and refine v sin U for WASP-10. We note that a larger stellar 
radius induces a greater planet size, which, in turn, could help to 
reduce the mass of WASP-10 b's core in the internal structure 
models. 

The durations of the individual transits were found to be the 
same for all four runs that shows that the ingresses or egresses 
of light curves are unaffected by occultations of starspots located 
close to the stellar disk limb. Limiting the O-C diagram to points 



based on high-quality photometry, it is obvious that the linear fit 
is an unsatisfactory solution because some data points deviate 
by up to 8cr. We re-analysed high-quality light curves available 
in the literature, i.e. from epochs 99 (Johnson et al. 20l0j) and 
242 (Maciejewski et al. 1201 II) . to measure the variation in the 
transit duration d tr . Following the rule of selecting the largest er- 
rors, the MC errors were taken as the final ones for the epoch 99 
light curve. In the second data set, a linear trend, clearly visible 
during the out-of-transit phase, was fitted simultaneously. Here 
the prayer-bead errors were found to be greater and taken in fur- 
ther calculations. We found d tr equal to 115.9^2 and 116.0^^ 
min for epoch 99 and 242 data sets, respectively. Both values are 
close to our determinations, a result that allows us to exclude an 
ingress/egress deformation by an occulting starspot as a source 
of the TTV signal. We also note that the redetermined mid-transit 
times were found to be consistent with the literature values and 
to lie well within the lcr error bars. 



5. Conclusions 

The signatures of WASP-10's activity have been observed in our 
high-precision transit light curves. Our analysis of system pa- 
rameters, which takes stellar activity into account, has indicated 
that the radius of WASP-10 b is smaller than most prior stud- 
ies report. This planet is three times more massive than Jupiter 
but has a similar radius. Its heavy-element mass, which is be- 
tween 270 and 450 M e , needs to be considered in the internal 
structure models to reproduce the planetary radius. Our new data 
points confirm that transit timing cannot be explained by a con- 
stant period if all literature data points are considered. They are 
consistent with the ephemeris assuming a periodic variation in 
the transit timing. We have found no means of explaining the 
observed TTV signal with light curve deformations caused by 
starspots on the stellar surface. This finding supports a scenario 
in which the second planet perturbs the orbital motion of WASP- 
10 b. However, we emphasis that additional high-quality transit 
light curves will be required to help us establish the true nature 
of this intriguing system. 
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